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Abstract 

The observed mass pattern of scalar resonances below 1 GeV gives preference to the 
tetraquark assignment over the conventional qq assignment for these states. We present a 
Oh| search for tetraquarks with isospins 0, 1/2, 1 in lattice QCD, where the isospin channels 1/2 
■ and 1 have not been studied before. We determine three energy levels in each isospin chan- 
nel using the variational method. The scattering states and possible tetraquark states are 
distinguished by considering the volume-dependence of spectral weights and by considering 
^ ! the time-dependence of correlators near t ^ T/2. We find no indication for light tetraquarks 
V/-^ \ at our range of pion masses 344 — 576 MeV. 
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> : 1 Introduction 

V . 

^ ■ In this paper we present a lattice study of the light scalar mesons cr, k, ao(980) and /o(980), 
^ ■ whose interpretation in terms of qq or tetraquarks qqqq is still not settled. 

The resonances ao(980) and /o(980) are well established experimentally [1]. The exper- 
imental evidence for a (or /o(600)) has also become strong in recent years [H [21 [3]. The 
position of the a pole has been determined from the experimental data using the model 
independent analysis with Roy equations, leading to mg. = 441 MeV and T„ = 54:4^^^^ 
MeV [3]. The position of the k (or ii'Q(800)) pole has been determined using a similar anal- 
ysis, leading to = 658 ± 13 MeV and = 557 ± 24 MeV [1]. Indications for k are also 
obtained from D, J/\I' and r decays [HIS], but the existence of k has still not been confirmed 
beyond doubt [1]. Assuming that k, exists, these states could form a SU{3) flavor nonet 
below 1 GeV illustrated in Fig. [21 In addition, there is a experimentally well established 
scalar nonet above 1 GeV as well as a few scalar glueball candidates. 
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The problem of identifying all observed scalar mesons as conventional qq states has a long 
history. The proposal that the states below 1 GeV could actually be tetraquarks qqqq dates 
back to Jaffe's work in 1977 [B] and this possibility is becoming ever more vivid in recent 
years: there has been a large amount of work based on phenomenological models [21 [6], [71 El [9] 
with reviews in [10]. In this paper we explore the possibility that the nonet below 1 GeV 
corresponds to Jaffe-type tetraquarks with a diquark anti-diquark structure [ggjs^^g^ [^^]3y,3c- 
We have the scalar (or "good") diquark [g^jg^^s^ in mind, which is anti-triplet in flavor and 
color and which is the most tightly bound among possible diquarks [TDl [TT]. Combining a 
scalar diquark with a scalar anti-diquark [qq\zj;i^, one gets a flavor nonet (3/x3/ = l/-|-8/) 
of color singlet (3c x 3c — ^ Ic) scalar particles. The nonet is expected to be light since the 
(anti) diquark is tightly bound and since this is a L = state (unlike qq which requires 
L = 1). The flavor pattern and the expected mass ordering for such a tetraquark nonet 
is schematically shown in Fig. (H along with the conventional qq nonet and the observed 
spectrum. 
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Figure 1: Schematic illustration of the observed spectrum for scalar mesons below 1 GeV (left), 
together with the expected mass spectrum for the nonet of scalar tetraquarks (middle), compared 
with a typical qq spectrum (right). 

Let us list just a few arguments which support a tetraquark interpretation of states below 
1 GeV, and disfavor a qq interpretation: 

• Most quark models place the ground state scalar qq with L = 1 above 1 GeV, along 
with well established axial and tensor qq mesons with L = 1. 

• The observed mass ordering < mao(980) can not be reconciled with conventional 
us and ud states since rrius > niud is expected due to > m^. This is probably 
the key evidence pointing to the tetraquark interpretation, where the J = 1/2 state 
has the flavor structure [ud] [ds] while the 1 = 1 state has the flavor structure [us] [ds] . 
The valence ss pair in a / = 1 tetraquark naturally makes this state heavier than the 
1 = 1/2 state and the resemblance with the observed spectrum (m^ < mao(980)) speaks 
for itself. 
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Although scalar mesons above 1 GeV are not the focus of this paper, let us mention 
that the observed ordering mi^o(i430) < ''^ao(i450) is also hard to reconcile with pure qq 
states. It has been argued that these states are mainly qq with a small admixture of a 
tetraquark state (via t'Hooft vertex), which leads to the observed mass ordering and 
decay pattern [8]. 

• The ao(980) is experimentally known to couple well with KK, since this decay channel 
is sizable in spite of the fact that most of the phase space for —>■ KK is closed. This 
again gives preference to the tetraquark over the qq interpretation: [Ms][(is] — > KK 
can proceed via "cheap" quark rearrangement [7] while ud KK is suppressed by 
the creation of ss pair. 

• Recently the interesting resonances X, Y and Z containing charm quarks were discov- 
ered and some of those are serious candidates for tetraquarks. If the BELLE discovery 
[12] of charged Z+(4430) in the channel ip'n'^ is confirmed, one has a clear indication 
for a four-quark state ccud. There have only been a few lattice studies using tetraquark 
interpolators aimed at these states [TB] . 

We note that there are also other proposals for solving the scalar meson puzzle, for example 
that light iso-scalars are largely glueballs [2] or that ao(980) and /o(980) are KK molecules, 
but these interpretations don't naturally lead to a nonet below 1 GeV. 

After numerous indications from experiments and from phenomenological models in favor 
of light scalar mesons as tetraquarks, it is very desirable to look for indications of tetraquarks 
in lattice QCD. The most straight forward quantities for a lattice study are the energies of the 
states with J^'^ = 0^~^ and with quark structure qqqq of a given flavor. The energy levels 
are extracted from the corresponding lattice correlation functions. The extracted energy 
level corresponds to the mass of a tetraquark state, if this exists. But there also unavoidably 
appear the energy levels corresponding to other 0"'"+ states of the same flavor, for example 
TTTT scattering states in case of a tetraquark with 7 = 0. In order to establish 

the existence of a tetraquark on the lattice, one has to identify all the relevant scattering 
states and distinguish these scattering states from a tetraquark state. We determine three 
energy levels for each isospin channel using the powerful variational method P^[TB| [T7]. The 
tetraquark and scattering states are distinguished by considering the volume- dependence of 
spectral weights and the time-dependence of correlators near t ~ T/2. 

There have only been a few lattice simulation of light scalar tetraquarks [HI [191 1201 l2T] . 
which are briefly reviewed in Section |H All of them consider only the 7 = state, while we 
consider the full flavor nonet with 7 = 0, 1/2, 1. All previous simulations except [IH] consider 
only the ground state, which is most likely the vrvr scattering state. In [19] the two excited 
states are extracted from a single correlator using the sequential empirical Bayes method [22] 
and the first excited state appears to be (t(600) with a tetraquark structure. This interesting 
result needs a confirmation using a different method, for example the variational method 
used here. 

Our main result is presented in Fig. [2l where the lowest three energy levels in each 
isospin channel are shown. We interpret the ground states in all channels as scattering states. 
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while the two excited states are above 2 GeV and can not correspond to hght tetraquarks. 
Therefore we find no indication for hght tetraquarks at our m„ ^ corresponding to m^r = 
344 — 576 MeV. However this does not exclude the possibility of finding tetraquarks in 
future lattice simulations with lighter quark masses or a more suitable basis of interpolators. 

Details of our simulation are presented in the next Section, while Section [3] provides our 
results. These are compared with previous lattice results in Section HI The last section 
provides conclusions. 
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Figure 2: The symbols present the three lowest energy levels from tetraquark correlators in 
/ = 0, 1/2, 1 channels at lattice volume 16^ x 32. The lines give analytic energy levels for scattering 
states: full lines present non- interacting energies ([7]), while dashed lines take into account tree-level 
energy shifts Appendix B explains why we don't see the energy corresponding to the scattering 
state Pi(l)P2(-l). 



2 Lattice simulation 

2.1 Correlation matrix and interpolators 

We calculate tetraquark correlation functions on the lattice, where a state qqqq with 0''""^ is 
created at t = and annihilated at some later time t. In each isospin channel / = 0, 1/2, 1 we 
use interpolators with correct quantum numbers and three different smearings i,j = 1,2,3 
that will be defined in ([5]) below 

C^.(t) = (0|O/(t)Of (0)10^=0^ = X](0|O/(f,t)Of (0,0)|0) . (1) 

X 
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The correlation matrix will be used to extract energy levels of the tetraquark system with 
the total momentum p equal to zero. 

We use theoretically well motivated diquark anti-diquark interpolators [SI El El [ID] 
composed of a scalar diquark in 3c 

[qQ]a = eabclqbC-f^Qc - QjC'y^qc] (2) 

and scalar anti-diquark in 3c 

[qQ]a = (!-abc[qhCi^Ql - QbC-f^qJ] , (3) 

where a, b, c are color indices. The [ud], [us], [ds] have flavor transformations just like s, d, u, 
respectively, while [ud], [us], [ds] transform like s,d,u. The combined tetraquark states 
[g(5][^'Q'] form a SU{3) flavor nonet and are expected to be light. For concreteness we 
simulate particular flavor combinations 

0^=0 = [ud][ud] , C^=^/2 _ ^ud][ds] , 0^=^ = [us][ds] (4) 

and assume = = m/. 

We use three different smearings of the quarks at the source Oj{0) and at the sink Oi{t) 
in order to extract more information from the system by means of the variational method. 
The three different tetraquark interpolators 0i=i,2,3 are constructed from smeared quarks of 
two different widths 

Oi = [qnQn][q'nQ'n] , O', = [qMlCQl] , O', = [qMWM (5) 

where g„ and refers to "narrow" and "wide" smearing, respectively. Smeared quarks 
are obtained from point quarks using spatially symmetric Jacobi smearing [23]. We use 
exactly the same two smearings (and same quark propagators) as applied in [21], which have 
approximately Gaussian shape and a width of a few lattice spacings (see Fig 1 of [21]). The 
same smearing parameters are used on different lattice volumes. 

2.2 Simulation details 

We use quenched gauge configurations, where the gauge fields are generated with the Liischer- 
Weisz gauge action [25] and are periodic in all four directions. All previous tetraquark 
simulations [THl [13, [201 [21] were quenched as well. As pointed out by Alford and Jaffe [TB] 
there is a good excuse to use the quenched approximation in this case: one is interested to 
look for a state with a well defined four quark content in this pioneering era of tetraquark 
searches on the lattice. The quenched approximation allows for this possibility as it does not 
mix qq <-> qqqq vac via sea-quark loops. We use another approximation, which discards 
this mixing: in the calculation of the correlator ([T]) we omit the single (b) and double (c) 
annihilation diagrams among Wick contractions in Fig. [31 

We work on two volumes V = L"^ xT = 16^ x 32 and 12^ x 24 at the same lattice spacing 
a = 0.148 fm, which is determined from the Sommer parameter ro — 0.5 fm. Using different 
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(a) (b) (c) 

Figure 3: Quark contractions for tetraquark correlators. 

volumes is crucial as it allows the determination of volume dependence for spectral weights. 
The quark propagators are computed from the Chirally Improved Dirac operator [26] with 
periodic boundary conditions in space and anti-periodic boundary conditions in time. We 
use mia = mu^a = 0.02, 0.04, 0.06 corresponding to = 344 — 576 MeV. The strange quark 
mass is fixed to m^a = 0.08, which leads to a vector meson mass closest to the physical 
(see Fig. 5 of [21] )• The resulting pion and kaon masses are given in Table [H The analysis 
is based on 96 configuration at 16'^ x 32 and 100 configurations at 12'^ x 24. 

2.3 The variational method and the expected eigenstates 

Next we discuss the extraction of the energy levels from the 3x3 correlation matrix ([T]). It 
decomposes in terms of n = 1, ..A^ physical states |n) as follows 

C^,{t) = Y.{0\O^An){n\O'/\0) e-^^' (6) 

n 

and the large-time correlation functions are dominated by the ground state. One of the 
physical states |?7,) is a tetraquark if it exists. In addition there are scattering states that 
unavoidably couple to our interpolators (jlj): vrvr for / = 0, Kir for 1 = 1/2 and KK, nrjgs 
for 1 = 1 {riss denotes a pseudoscalar meson composed of ss). On a lattice with spatial 
extent L and periodic boundary conditions for (anti) quarks, the energies of scattering states 

— * — * 

Pi{k)P2{—k) with total momentum zero come in a tower 

^AO)P.(-.) = + ... , + ^m^,+ (^y, ... (7) 

for integer j in assuming non-interacting pseudoscalars. The lowest few energy levels are 
well separated at the volumes we are using. In the interacting theory, the energy levels in a 
finite box are shifted by Liischer shifts [27] . These are expressed in terms of infinite- volume 
scattering lenghts and are given by [28] 

dE'::U = -Vi^flL') dE^,^,,, = -l/{flL^) (8) 

at tree-leve0. 

^Liischer [57] provides sliifts in full QCD. The tree-level shifts in quenched QCD agree with full QCD. 
The quenched shifts beyond tree-level have only been derived for the case of tttt scattering [29] . 
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In order to show the existence or absence of tetraquarks on the lattice, one has to estabhsh 
energy-levels for a tetraquark as well as for the few lowest scattering states. Since fitting the 
sub-leading exponentials in is very unstable, we use the variational method. We therefore 
compute the eigenvalues and eigenvectors which diagonalize the 3x3 hermitian correlation 
matrixl^ 

C{t)v^it) = Xn{t)Vn{t) (9) 

at each time slice. The eigenvalues are dominated by a single physical state [12] 

A„(t) = ^„e-^"* (1 + 0(e-^"*)) (10) 

which can be easily understood if there were only = 3 physical states. The corrections 
are given by the energy difference to the nearest energy level A„. The eigenvectors Vn(t) are 
orthogonal and represent the components of physical states in terms of the chosen three- 
dimensional variational basis (Q. We also have determined the energy levels using the 
generalized eigenvalue problem C {to)~^^'^C {t)C {to)~^^'^v!^{t) = X'j^{t)v'^{t) for various values 
of to and the extracted energies agree with the results from the standard eigenvalue problem 
([9]) withing the errors. In the following we only present results based on the standard 
eigenvalue problem, where the determination of w as a function of L is more direcllf|. 

All eigenvalues are fitted using uncorrelated fits, which explains a rather small /d.o.f in 
the Tables. All quoted errors are statistical and are determined using the jackknife method. 



3 Results 
3.1 1=0 

Energies 

In this channel we need to identify the vrvr scattering states and find out if there is an 
additional candidate for a with the content [ud] [ud] . 



^ All our correlators Cy (t) are symmetric with respect tot T—t and we fold them as i [Cy [t)+Cij [T—t)] 
before computing the eigenvalues. 

•^In the generalized problem C(to) defines the normalization and A'j(io) = 1 by construction. 
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three eigenvalues for 1=1, m^= 0.04, V=16^x32 
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Figure 4: Effective masses for the three eigenvalues Ao,i,2(i) in / = 0, 1/2, 1 channels for mia = 0.04 
and V = IQ^ x 32. The full diamonds were obtained using rrieff defined in (jl9|2ip . while others 
were obtained using conventional cosh-like definition (jl4p . The lines give energy levels for scattering 
states: full lines present non- interacting energies ([7]), while dashed lines take into account tree-level 
energy shifts ([8]). g 
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Figure 5: Eigenvector components Re{vl^) as a function of t for three eigenstates n = 0,1,2 in 
terms of the basis Oi ([5]). We find Im{v\^) ~ 0. The figure shows the case 1 = 0, m/o = 0.04 
and V = 16^ x 32 and others are very similar: eigenvectors show almost no isospin, quark mass or 
volume dependence. 

The three lowest energy levels are given by the effective masses for the three eigenvalues 
\n{t) at large t. They are plotted in Fig. for the case of mi = 0.04 and V = 16^ x 32. 
The corresponding eigenvectors, given in Fig. [5l illustrate the decomposition of the three 
eigenstates in terms of our basis ([5]): \n) = The energy of the ground state is 

close to 2m^, indicating that the ground state is a 7r(0)7r(0) scattering state. Below we will 
present two other arguments which support this interpretation of the ground state. The 
most important finding is that there is a large gap above the ground state: the first and the 
second excited states appear only at energies above 2 GeV. Whatever the nature of these 
two excited states is, they are much to heavy to correspond to cr(600), which is the light 
tetraquark candidate in this channel. The two excited states may correspond to ir{k)7r{—k) 
with higher k or to some other energetic state. We refrain from identifying them with certain 
physical objects since such energetic states are not the focus of our present study. 

A very similar conclusion applies for the other light quark masses and for the other 
volume. All the ground state energy levels are close to 2m^, while the excited energy levels 
i?i 2 are above 2 GeV, as shown in Fig. [2] and Table [2] 0. Apart from 7r(0)7r(0), we find no 
state with energy close to m„, so we find no indication for a light tetraquark in the 1 = 
channel. Since all our pion masses are just above 300 MeV, this is not in conflict with the 
simulation of the Kentucky group [19], which finds indication for a tetraquark with mass 
~ 550 MeV at = 180 - 300 MeV (but not above that). 

Looking at the spectra in Figs. H] and [2] the question araises why there is no state close 

■^The fitting forms for the extraction of energies are presented in the next subsection and in Appendix A. 
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to the energy 2^ml + (27r/L)2 of a -7r(l)7r(— 1) state. In Appendix B we argue that the 
tower of few lowest scattering states 7i{k)n{—k) contributes to the lowest eigenvalue, which 
explains why its effective mass is not flat at intermediate t. This is due to the fact that 
our basis (jS]) does not disentangle the few lowest scattering states into separate eigenvalues. 
They would contribute to separate eigenvalues by using a larger basis or using interpolators 
with a definite momentum projection for each pion. 

Time-dependence of the eigenvalues 

We were surprised to find that the ground state eigenvalues do not have a conventional 
time-dependence iy[e~^* + e"^'-"'""*-*] near t ~ T/2. This was first noticed by looking at the 
cosh-type effective mass (1141) . which is decreasing near t ~ T/2 (see empty symbols in Figs. 
|4]and[8]). In Appendix A we derive that the eigenvalues for the vrvr state receive a constant 
in addition to the conventional term w[e~^^ + e"^'-"^"*-'] (fT8|) . Our data for X^'^'^^t) agrees 
well with this analytic expectation and the effective mass (fT9l) which takes this into account 
has a plateau near t ~ T/2 (see full symbols in Figs. HlandlH]). The non-conventional time- 
dependence of the ground state is the second indication that this is a scattering state. The 
values of Eq and wq in Table [2] are obtained from the three-parameter fit of Ao(t) f|T8l) . The 
presence of an additional constant increases the uncertainty of the resulting Eo and this is 
the main reason that we refrain from studying the nir energy shifts. 

All excited states are obtained from the conventional fit to Ai 2(t) = u'i^2[e~^^'^* + 

g-Sl,2(T-t)j_ 

Spectral weights 

Let us describe a method, which considers w{L) in order to allow the distinction between 
the one-particle (tetraquark) and the two-particle (scattering) states on the lattice. 

The L dependence of the spectral weight w for a scattering state P1P2 has been derived in 
coordinate space in Section III E of [SO]- We re-derive it in momentum space by evaluating 
the contribution in Fig. [6^ for finite L, infinite T and p = 



{OMPiP2){PiP2\Ol{0)) oc / dp^ e 



dk dh \{0^\Plik)P2i-k))\^ 
(27r)3 27 [kj + P + mj] [{p^ - k^f + P + ml 



j 12 

The leading volume dependence of the coefficient in front of the exponent (called the spectral 
weight w) for a given physical state j is 



and it comes from dkj = 271 /L. The energies E^^^ = \J m? + (27rj'/L)2 in flTT]) are only 
mildly dependent on L. The couplings {Oi\Pi{j)P2{~j)) at given j are volume independent 
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when the size and the shape of interpolators do not depend on L (which is true for our 
interpolators) and when the coupling is non-derivative. In the case of derivative coupling 

[8], which is proportional to kpjip^^ = \J m\ + {2T{j/LY\Jm\ + {2nj/Ly + {27rj/Ly, there 

is a mild dependence on L for j ^ 0, but this dependence is sub-leading with respect to the 
dependence (fT2!) . 

The contribution of a one-particle (tetraquark) state in Fig. [6)d gives 

{OMT){T\Olm oc / dp, e^^' ^— ^ (X (13) 

and the spectral weight is expected to be almost independent of the volume. 

P2(p-k) 




(a) Pi(k) (b) 



Figure 6: Diagrams relevant for deriving the volume dependence of spectral weight for two-particle 
(a) and one-particle (b) state contributions. 
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Figure 7: The ratio of spectral weights wq{L = 12)/wq{L = 16) for / = 0, 1/2 as computed from 
the ground state eigenvalues for two volumes L^. 

Our spectral weights w are presented for both volumes in Table [2l The ratio Wq{L = 
12)/wq{L = 16) for the ground state is plotted in Fig. [71 where the individual errors 
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are summed in quadrature A{wi2/wiq) = ^/\Awi2/wiq\'^ + \Awiq Wi2/wIq\'^. The ratio is 
consistent with 16^/12^ within large error bars and its w oc dependence is the third 

indication that our ground state is a vrvr scattering state. 

3.2 7 = 1/2 

The conclusions regarding 1 = 1/2 are very similar to / = 0. 

The ground state energy in is close to m7r + mx (see Figs. Il]and[2]), which indicates it is a 
7r(0)Ji (0) state. The corresponding eigenvalue has a non-conventional time dependence near 
t ^ r/2, which agrees with the analytic expectation (!20|) for a 7r{0)K{0) state. The cosh-type 
effective mass does not have a plateau near t ~ T/2, while the corrected effective mass (12T!) 
has a plateau (see Fig. [S]). The three parameter fit of \o(t) (12U]) gives Eq and Wq in Table [31 
The volume dependence of the spectral weight for the ground state Wo{L = 12) /wq{L = 16) 
is consistent with 16'^/12^ (see Fig. [7]) which is another indication for a Ktt scattering state. 

The energies of the first and the second excited states in Figs. H] and [2] are above 2 GeV. 
So we do not find an indication of a tetraquark with mass close to the mass of k(800). 

3.3 1 = 1 

There are two scattering states KK and irriss with similar energies at / = 1, which makes this 
channel more challenging. The major features of the results are the same as for J = 0, 1/2: 
the ground state energy in Figs. [2] and H] is close to 2mK and -|- m^^^ = -|- ^1^^^^ , while 
the excited states are above 2 GeV. Again, we do not find an indication for a tetraquark 
with mass close to the mass of ao(980). 

However, in this channel our conclusion is not so firm since we were not able to analyze 
Ao(t) using a two-state fit with KK f|T8l) and vrr/s^ fl20l) . Such a fit is not stable with our 
data, especially since the two states are close in energy and since there are additional terms 
due to finite T in fll8ll20p . We present only approximate ranged for the ground state energies 
in Table H] and Fig. [2], while we are unable to present reliable values for the corresponding 
spectral weights. 

Due to the presence of two towers of scattering states in addition to a possible tetraquark, 
the 7 = 1 channel remains a challenge with the methods available at present, even for future 
simulations. 

4 Previous lattice simulations of light tetraquarks 

Finally we compare our results with the results of previous simulations [TSl dSl [201 [2T] . Like 
ours, all these simulations are quenched and omit the contributions of annihilation diagrams 

^The naive cosh- type effective mass ([T4)) for V ~ 16"^ x 32 has a brief plateau at t = 8 — 10 and a sizable 
fall-off for t = 10 — 16. We give approximate ranges of ground state energies as rn^-^^ — ?^^g^J in Tableland 
Fig. [21 which should serve only as guidance. 
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(b), (c) in FigO In contrast to the present work, they all consider only J = |^ and only 
the ground state (except for [T^). 

Alford and Jaffe [H] computed the ground state energy Eo for tttt interpolators at fixed 

~ 800 MeV and various L. They find that the energy shift EQ — 2m-,^ does not completely 
agree with the (full ChPT) analytic prediction for tttt scattering at finite volume. They 
interpret this as a possible indication for a tetraquark. 

The Kentucky group [19] computed the tttt correlator using overlap quarks for = 
182 — 760 MeV and two volumes L^. Using the sequential empirical Bayes method [22] 
they were able to extract three states from a single correlator. The ground state with 
Eq ~ 2m^ + dE{L) and wo{L = 12)/wo{L = 16) ~ 16712^ is interpreted as 7r(0)7r(0). The 
first excited state with Ei ~ 600 MeV and wi{L = 12)/wi{L = 16) ^ 1 gives indication for a 
with a tetraquark structure. They find this state only for m,r < 300 MeV, so their results are 
not in contradiction with our results, which give no evidence for a tetraquark at > 300 
MeV. Their second excited state has E2 — E^^^^^^~^^ and is interpreted as 7r(l)7r(— 1). This 
is one of the few available lattice studies which was able to separate the one and two-particle 
contributions and gives impressive evidence for a as a tetraquark. However, it relies on the 
ability of the sequential method to extract three states from a single correlator and therefore 
needs confirmation using other methods. 

The Japanese group [2U] determined the ground state energy for a diquark-antidiquark 
interpolator [M(i]['U(i] at mu^d = [nis, 2ms]. They find that the extracted energy is close to the 
TTTT energy in the case of periodic as well as in the case of hybrid boundary conditions. They 
conclude the ground state is not a tetraquark state. 

The authors of [21] use the variational method and extract only the energy of the ground 
state, which is close to 2m^ for their range of pion masses. 

5 Conclusions and outlook 

The question whether the observed scalar resonances ao(980), /o(980), cr(600) and k(800) 
correspond to conventional qq or tetraquark nonet is still not settled. The observed mass 
pattern gives preference to the later. 

We presented a search for possible light scalar tetraquarks by means of a quenched 
lattice simulation. In each isospin channel we extracted the three lowest states from a 3 x 3 
correlation matrix with diquark-antidiquark interpolators. The energy levels are shown in 
Fig. [2], which is the main result of our paper. The ground states are found to be consistent 
with the scattering states. The first and second excited states have energies above 2 GeV, 
so they can not correspond to resonances below 1 GeV. 

In conclusion, we find no indication for light tetraquarks at our range of pion masses 
344 — 576 MeV. However, one should not give up hopes for finding these interesting objects 
on the lattice. Indeed, our simulation does not exclude the possibility of finding tetraquarks 

^The authors of [T8l[2T] also consider the 1 = 2 channel, but this does not belong to the nonet considered 
here. 
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for lighter rriu^d or for a larger (different) interpolator basis. A stimulating lattice indication 
for cr as a tetraquark state at = 182 — 300 MeV has already been presented in [T^ . 

The present and past pioneering quenched tetraquark simulations, which discard annihi- 
lation diagrams, provide valuable information on the states with a definite quark assignment. 
The final conclusions will have to await dynamical simulations incorporating both annihila- 
tion quark diagrams and the qqqq ^ qq ^ vac mixing. 
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A Time dependence of correlators for scattering states 
at finite T 

In this Appendix we discuss an interesting observation, related to scattering states on the 
lattice with finite time extent T. We noticed an unconventional behavior by looking at the 
cosh-type effective mass 

(14) 



for the ground state, which keeps falling even at large t close to T/2 (see empty symbols in 
Figs. |l]and[8]). This indicates that Xo(t) does not have the conventional time-dependence 

We believe this is due to the fact that the ground state is a scattering state PiP2- Let us 
derive the time dependence of a correlator C{t) which creates and annihilates a state by the 
interpolator O. We consider the case where the lowest physical state coupling to O is the 
scattering state PiP2- An example is P1P2 = 7i~^K^ with O = [ud][ds] ot O = 7r~^K^, which 
is relevant to our correlators with 7 = 1/2. We start from the basic definition of C(t) [31] 



at) = lTr[e-^^C)(t)C»t(0)] =|Tr[e-^(^-*)Oe-^*C»t] 
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When t and T — t are large, only the states n, m = 0, P1P2, Pi, P2 contribute and the non- 
vanishing terms in this limit are 

C{t) = ^(^(0|e-^(^"*)C|PiP2)(Ai'2|e~''*C^|0) + (PM|e-^(^-*)O|0)(0|e-^*C^|PM) 
+ (P/|e"^(^-*)0|P2)(P2|e-^*Ot|pt) + (p||e-^(^-*)0|Pi)(Pi|e-^*Ot|pt)^ 



|(^|(0|O|PiP2)|V^"^"^ *+ |(P/P2^|O|0)|V^"^"^(^"*) 
+ |(P/|0|P2)|V^^i('^-*)e-^^2* + |(pj|0|Pi)|2e-^^2(^-*)e-^^i*^ 



(16) 



where H\0) = 0. In first term Pi and P2 propagate forward in time, in the second they 
both propagate backward, while in the third and fourth term one propagates forward and 
the other backward. For the case of our ground state Pi(0)P2(0) with the anti-periodic 
propagators in time this gives the appropriate form to fit the eigenvalue 

;^Pi(o)P.(o)^^^ = u;[e-^«* + e-^«(^-*)] + A[e-'""*e-'"^^(^-*) + e— ^^^e-'""^^-*)] (17) 

with E'o ~ mpi -|- mp2- A similar relation was provided without derivation in the Appendix 
of [32]. 

Therefore we fit the eigenvalue for the ground state 7r(0)7r(0) to 

Xl-{t) = i/;o[e~''°* + e-^°(^-*)] + const (18) 

at large t and the resulting fit parameters wq, Eq and const are given in Table [21 All 
terms in (ITTl) are equally important near t ~ T/2 since (e~'^°* -|- e~^°^'^~^^)t=T /2 / const = 
2e~^'^'^/'^ / const is of the order one, as shown in Table [2l The effective mass computed after 
subtracting the resulting constant (|T8l) 

t + 1/2 , t + 1/2 irr, ,N 

Xl'^it)- const _ e-^^^ff * + 6"™=// 



Aj^(t + 1) - const ^-miy/^ (t+i) + e-</)'" (^-t-i) 

is plotted by full symbols in Figs. H] and [H It has a plateau at large t in contrast to 
the effective mass computed from f[T4l) . which is another indication that the ground state 
eigenvalue corresponds to vrvr. 

The 7 = 1/2 ground state 7i{0)K{0) is fitted to 

A-^(t) = wo[e-^°' + e-^o(^-*)] + ^[e-™-*e-'"^(^-*) + e—^'^e-'"-^^-*)] (20) 

and the resulting fit parameters Eq, wq and A are given in Table [31 Note that A is comparable 
to Wq, so all four terms are of similar size near t ^ T/2. The and mx were fixed to the 
measured values in (l20l) and we verified that the variation of the results is negligible if m-jr^K 
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are varied in the ranges given in Table [H The effective mass computed after subtracting the 
last term in (I2U]) 



A^^(t) - yi[e-'»-*e-'"^(^-*) + {t^T-t}] 
X^^it + 1) - A[e-™-(*+i)e-™^(^-*-i) + {t^T-t}] 



e 



(21) 

has a plateau at large t (see Figs. HI [S]), which is another indication that ground state 
eigenvalue corresponds to the ttK scattering. 
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Figure 8: Effective masses for the / = 0,1/2 ground states at various m/o and y = 16^ X 32. 
Empty symbols are obtained using "naive" cosh- like definition for m(,ff (jl4p . while full symbols 
take into account the correct form: ([T9j) for tttt scattering and (j2ip for i^Tr scattering. 
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B Variational analysis with several states contributing 
to a single eigenvalue 

A distinct feature of our spectrum is that there is no eigenstate corresponding to the first 
excited scattering state Pi{l)P2{—l) (see Fig. Hj). Another feature is that the effective mass 
for the ground state is not flat at intermediate t and seems to contain important contributions 
from several physical states, not just Pi(0)P2(0) with a correction flTU]) . These observations 
indicate that our interpolator basis (Q does not decouple the few lowest scattering states to 
separate eigenvalues and several of them contribute to the ground state eigenvalue. 

Let us present an example for a specific mechanism that may be responsible for that. 
The point source couples to all the scattering states equally 

[QQ][QQ]point = Cpoint^9k\Pi{k)P2{-k)) + ... , (22) 

k 

up to a factor gj:, which gives a Lorentz structure of the coupling ((7^ = 1 for non- derivative 
coupling and (7^ = ki^k2 for derivative coupling [8]). For a spatially extended source one 
expects additional fc-dependence foi{k) in (122!) given by the shape of smearing. Our inter- 
polators ([5]) are not point, but are all rather narrow with an extent of a few lattice spacings. 
Let us explore the consequences of the approximation that all our sources behave close to 
point-like. In this case the few lowest scattering states couple to a given interpolator equally, 
but there is a different overall magnitude for each interpolator 

n \Pl{k)P2{-k)) + annnn\a) + KnnnW) 

k 

k 



k 



Here we assumed that our interpolators couple only to two additional physical states a and 
h. Given these linear combinations, one can construct the corresponding 3x3 correlation 
matrix and it can be easily shown that its eigenvalues are 



= ^0 ^9k 
K{t) = Wae 



g_(£;Pi(fe)+£;P2W)t 

i Wo 



k 

-Eat 



^Pl{k)^P2{k) 



mp,mp2 E^^WEP^m + - 



Xbit) = Wb e-^'* . (23) 

The physical states a and b get their own exponentially falling eigenvalues, while a whole 
tower of scattering states contributes to a single eigenvalue in this approximation. Correc- 
tions to (1231) come from additional heavy physical states and from slight /c-dependence of 
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shape functions foX^) for our sources. We note that a related variational problem has been 
studied in [33], where states with non-exponential time dependence (ghosts) were found to 
contribute to a separate eigenvalue. 

ground state for 1=0, V= 16^x32 

2 
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0.6 
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t 
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t 

Figure 9: Symbols present effective masses (jl9p for the ground state eigenvalues in / = 0,1/2 
channels. Lines are effective masses for analytic prediction (j23p where a tower of tttt (/ = 0) or ttK 
(/ = 1/2) scattering states contributes to Xo(t). We sum (p3|) over five lowest and verify that 
the contribution of higher is negligible in the plotted time range. We use gj: = 1, the measured 
values of TTiTT^A- and the ground state energies mp^ +mp^ (full lines) or mp-^ +mp^ + dE*'^^^ (dashed 
lines) . 

We find that the described scenario describes our ground state eigenvalues rather well. 
Fig. [9] shows that the effective masses ( fT9l[2Tl) agree quite well with the effective masses 
for the analytic prediction Ao(t) in fl23|l with g^: = 1. We conclude that (i) a coupling 
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[qq\ [qq\ P1P2 is close to a non-derivative coupling with = ^ for our interpolators and 
(ii) a tower of few lowest scattering states vrvr, Kn is indeed contributing to our ground state 
eigenvalues for / = 0, 1/2. 
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mia 


niga 




[MeV] 


rriKCi 


rriK [MeV] 


fit-range 


0.02 


0.08 


0.259(3) 


344(4) 


0.397(3) 


528(4) 


7-12 


0.04 


0.08 


0.357(3) 


475(4) 


0.433(3) 


576(4) 


7-12 


0.06 


0.08 


0.433(3) 


576(4) 


0.466(3) 


620(4) 


7-12 


0.08 


0.08 


0.500(3) 


665(4) 


0.500(3) 


665(4) 


7-12 



Table 1: The pion and kaon masses from our quenched simulation at V = 16^ x 32 and = 1.33 
GeV, obtained from the ground state eigenvalue of 3 x 3 correlation matrix presented in [24]. The 
values at 12^ x 24 agree within the errors. 



/ = 


mia 


V 




Wn 




fit-range 




const 


2e-^oT/2 




d.o.f. 


const 




0.02 


16^32 


0.54(3) 


1.8(6) ■ 


1015 


9-16 


0.01 


6.1(16) 


10" 


1.0 


ground 


0.04 


16^32 


0.71(2) 


5.1(13) 




11-16 


0.002 


5.8(13) 


. 10^ 


2.0 


state 


0.04 


12^24 


0.72(5) 


1.3(5) ■ 


10^5 


8-12 


0.04 


2.6(11) 




1.7 




0.06 


16^32 


0.88(2) 


4.4(7) ■ 


10^4 


11-16 


0.02 


3.1(4) 


10^ 


2.3 




0.06 


12^24 


0.86(3) 


8.8(23) 


10^4 


8-12 


0.03 


1.9(9) . 


lO^o 


3.0 




0.04 


163 32 


1.63(5) 


1.1(2). 


1013 


3-8 


0.07 








first 


0.04 


123 24 


1.57(9) 


1.1(3)- 


10^3 


3-7 


0.07 








exc. 


0.06 


163 32 


1.70(4) 


9.9(13) 


.10^2 


3-8 


1.5 








state 


0.06 


123 24 


1.69(6) 


1.1(2). 


10^3 


3-7 


0.1 










0.04 


163 32 


1.88(5) 


4.4(7) . 


10^1 


3-6 


0.1 








sec. 


0.04 


123 24 


1.96(8) 


5.6(14) 


.10^1 


3-6 


0.4 








exc. 


0.06 


163 32 


1.97(4) 


4.4(5) . 


10^1 


3-7 


0.2 








state 


0.06 


123 24 


1.97(5) 


4.5(8) . 


10^1 


3-7 


0.7 









Table 2: The energy levels En and spectral weights Wn from three eigenvalues in / = channel, for 
various m/ and V. The energies in physical units are obtained by multiplying Ea with ~ 1.33 
GeV. The ground state eigenvalue is fitted to (jlSp . which contains additional const and its relative 
importance at t = T/2 is represented by 2e~^^/^/const. For the lowest mass mia = 0.02 some 
results are not presented as they are to noisy. 
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J = 1/2 


rriia 


V 


End 






fit-range 




A 




d.o.f. 




0.02 


16^32 


0.67(3) 


9.2(29) 


10^4 


10-16 


0.001 


5 Sfl 31 


■ lOi^ 


ground 


0.04 


16^32 


0.79(3) 


4.3(13) 


10^4 


12-16 


0.005 


2.0(5) ■ 


lOi^ 


state 
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12^24 


0.78(4) 
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X \J 
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sec. 


0.04 


123 24 


1.97(7) 
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lOii 
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4.5(7) ■ 
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0.7 







Table 3: Same as Table [2] but for 1 = 1/2 channel. The ground state eigenvalue is fitted to (j20p 
with an additional term which is proportional to the parameter A. 



I = 1 


mia 


V 


Ena 






fit-range 


x' 

d.o.f. 




0.02 


16332 


0.68-0.82 










ground 


0.04 


16332 


0.78 - 0.92 










state 


0.06 


16332 


0.88-1.00 
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1.65(4) 


1.1(2). 


1013 
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0.05 




0.02 


123 24 


1.57(10) 


1.1(3). 


1013 


3-6 


0.1 




0.04 


163 32 


1.69(4) 


1.0(1). 


1013 


3-8 


1.2 


first 


0.04 


123 24 


1.67(6) 


1.1(2). 
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3-7 


0.1 


exc. 


0.06 


163 32 


1.73(4) 


9.7(12) 
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3-8 


2.2 


state 


0.06 


123 24 


1.73(5) 
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1013 


3-7 


0.2 




0.02 


163 32 


1.94(5) 


4.6(8) . 


lOii 
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0.1 


sec. 


0.02 


123 24 


1.91(8) 


4.4(11) 
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0.4 
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0.04 


163 32 


1.97(4) 


4.4(5) . 


lOii 


3-6 


0.3 


state 


0.04 


123 24 


1.96(5) 


4.4(7) . 


lOii 


3-6 


0.8 




0.06 


163 32 


2.01(4) 


4.4(5) . 


lOii 


3-7 


0.3 




0.06 


123 24 


1.99(5) 


4.3(6) . 


lOii 
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0.9 



Table 4: Same as Table[2]but for 1=1 channel. We are unable to perform a two-state {KK, tttJss) 
fit for the ground state eigenvalue, so we provide only informative energy ranges given by the naive 
cosh-type effective mass rrv^-rr — m^'f • 
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